Photonic quasicrystals for general purpose nonlinear optical frequency conversion 
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We present a general method for the design of 2-dimensional nonhnear photonic quasicrystals 
that can be utilized for the simultaneous phase-matching of arbitrary optical frequency-conversion 
processes. The proposed scheme — based on the generalized dual-grid method that is used for con- 
structing tiling models of quasicrystals — gives complete design flexibility, removing any constraints 
imposed by previous approaches. As an example we demonstrate the design of a color fan — a non- 
linear photonic quasicrystal whose input is a single wave at frequency uj and whose output consists 
of the second, third, and fourth harmonics of u, each in a different spatial direction. 

PACS numbers: 42.65.Ky, 42.70.Mp, 61.44.Br, 42.79.Nv 
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The problem of phase matching in the interaction of 
light waves in nonlinear dielectrics became immediately 
evident as the first theories describing such interaction 
were developed 0. Put simply, nonlinear interaction 
is severely constrained in dispersive materials because 
the interacting photons must conserve their total energy 
and momentum. Even the slightest wave-vector mis- 
match appears as an oscillating phase that averages out 
the outgoing waves, hence the term "phase mismatch" . 
One approach for treating the problem uses the bire- 
fringent properties of specific materials and by playing 
with the polarizations of the interacting waves achieves 
phase matchingjdjs . A second approach, suggested over 
4 decades ago [l|, y and known today as "quasi-phase- 
matching" , is to modulate the sign of the relevant com- 
ponent(s) of the nonlinear dielectric tensor at the period 
of the oscillating mismatched phase thereby undoing the 
averaging. Quasi-phase-matching has been generalized 
from simple 1-dimensional periodic modulation to 2- 
dimensional periodic modulation d, S 0J3as well as 
1-dimensional quasiperiodic modulation Il2l , 
allowing greater flexibility in phase-matching multiple 
frequency-conversion processes within the same photonic 
crystal. Here we present the full generalization of the 
method that enables the design of nonlinear photonic 
crystals that can simultaneously phase-match any arbi- 
trary set of frequency-conversion processes in any spatial 
direction. This design flexibility is ideal for the realiza- 
tion of elaborate multi-step cascading effects 0, 0| , as 
demonstrated by the color fan example (Fig. ^ at the 
end of this article. 

To understand how the method works it is convenient 
to adopt the view taken in condensed matter systems. 
Recall that momentum conservation is a direct conse- 
quence of having continuous translation symmetry. In 
crystals, whether periodic or not, continuous translation 
symmetry is broken, and momentum conservation is re- 
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FIG. 1: The color fan. 



placed by the less-restrictive conservation law of crystal- 
momentum. The total momentum of any set of inter- 
acting particles in a crystal — whether they are electrons, 
phonons, or photons — need only be conserved to within 
a wave vector from the reciprocal lattice of the crystal, 
giving rise to so-called umklapp processes. Thus, all one 
needs to do is to construct an artificial photonic crystal 
whose nonlinear susceptibility is ordered on the appropri- 
ate length scale such that its reciprocal lattice contains all 
the necessary wave vectors, needed to phase-match any 
required frequency-conversion process. The linear sus- 
ceptibility should remain constant in space so that the 
dispersion relation will not vary from point to point. We 
show here a method — borrowed from the theory of qua- 
sicrystals ^3 — to construct ordered structures in real 
space with any given Fourier components. Such struc- 
tures are in general quasiperiodic, thus we are concerned 
here with the design of nonlinear photonic quasicrystals 
(for a definition of the term see the discussion in Ref.0). 

To be more specific, consider for example, two incom- 
ing waves (a;i,ki) and {uj2,'k.2) interacting in the plane 
of a 2-dimensional dielectric whose total area is A, via its 
lowest-order nonlinear susceptibility tensor x*-^-*, to pro- 
duce an outgoing wave (t^s, ks), with lo^ = cji -f a;2 and a 
wave vector mismatch Ak = ki -|-k2 — ka. Assuming neg- 
ligible depletion of the input waves, the amplitude of the 
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outgoing wave as it emerges from the interaction region 
A is 

EsiAk) =T [ g{Y)e'^^-^(fr, (1) 

J A 

where F is a constant that depends on the amplitudes of 
the incoming waves and on the indices of refraction of all 
three waves, and g(r) is the relevant component of x*-^', 
coupling the three waves. Clearly, if g is independent of r 
the oscillating mismatched phase averages out the output 
signal giving, for a rectangular A, a two-dimensional sine 
function that tends to a delta function (5(Ak) as the in- 
teraction area A tends to infinity. In general, the output 
signal i?3(Ak) is proportional to the Fourier transform 
of g(r)(y9A(r), where V'A(r) is the characteristic function 
of the interaction area A (equal to 1 if r G ^4, and oth- 
erwise). We can design this Fourier transform to peak 
at Ak by ensuring that Ak is a vector in the reciprocal 
lattice of g[r). 

Consider more generally a set of D nonlinear interac- 
tions of triplets of waves, coupled by x^^-* as in with D 
2-dimensional wave vectors Ak'^' [j = 1 . . . D], describ- 
ing their mismatches. We wish to design a nonlinear pho- 
tonic crystal ^(r) that will simultaneously phase-match 
all D interactions. If, by accident, all D mismatch vectors 
can be expressed as integral linear combinations of two of 
them, Ak*^^^ and Ak*^^\ one can view these two vectors 
as the basis of a 2-dimensional reciprocal lattice. One 
then finds the dual real-space basis-vectors a*^^^ and a*^^^ 
whose reciprocal vectors are Ak*^^^ and Ak*^^) through 
the condition 

a(^) . Ak(^'' =27r5y, (2) 

and generates a real-space lattice consisting of points at 
all integral linear combinations of a^^^ and a^^^. The 
Fourier transform of a set of delta functions at these lat- 
tice positions is a set of delta function at all the recipro- 
cal lattice positions generated by integral linear combina- 
tions of the two mismatch wave vectors, as required. This 
is the procedure that has been used until now |E S S ■ 
In general, when the mismatch vectors do not all be- 
long to a reciprocal lattice of a periodic crystal, we 
need to construct a quasicrystal whose Fourier trans- 
form contains all the mismatch vectors. To do so 
we use a generalized version of de Bruijn's dual grid 
method [191 used for the construction of tiling 

models of quasicrystals. We begin by selecting a set of D 
2-dimensional real-space vectors a'-'' [j — 1 . . . D], called 
tiling vectors. We use a generalization of the orthogonal- 
ity condition known as the Ho condition j2^ . 

D 

^a^AfcW =27r5^.. (3) 

Note that in the special case of a periodic structure the 
Ho condition (|2Jl is simply the completeness relation that 



follows from the orthogonality condition ^ used to de- 
fine the reciprocal vectors. If we were now simply to 
generate all points at integral linear combinations of the 
D tiling vectors we would get the unwanted outcome of a 
dense filling of real space. To avoid this situation we con- 
struct the dual grid, whose topology determines which of 
the integral linear combination of the tiling vectors are 
to be included in the tiling. 

To construct the dual grid we associate with each mis- 
match vector Ak^-') = 27rn'^-') / Lj an infinite family of 
parallel lines normal to the unit vector n^^^ , separated 
by a distance Lo , and shifted from the origin (in the di- 
rection of — n(^5) by an amount /jLj, for some chosen 
set of grid shifts < /j < 1. This grid is dual to the 
tiling in the sense that each intersection of lines in the 
grid corresponds to a tile in the tiling whose edges are 
the tiling vectors a*^-'^ associated with the families of the 
intersecting lines. In addition, each cell in the dual grid, 
labeled by D integers Uj , determines a vertex in the tiling 
at position njaS^\ where Uj is the number of lines of 
the j*'* family separating the cell from the origin. Only 
those linear combinations that correspond to cells in the 
dual grid are included in the tiling. 

The canonical choice of tiling vectors 0, section IV. C], 
which is convenient for the analytical calculation of the 
Fourier transform, is obtained by viewing the D 2- 
dimensional mismatch vectors as 2 D-dimensional vectors 
(Afc^^-* , . . . , A/c^^'') = 1,2] spanning a 2-dimensional 
subspace of a Z?-dimensional vector space. One can 
then choose D — 2 additional D-vectors (qjP, . . . , qjF^) 
[fi — 3...Z?], orthogonal to the first two, spanning 
the remaining {D — 2)-dimensional subspace. This pro- 
cess adds D — 2 components to each of the mismatch 
vectors extending them into the D D-dimensional vec- 
tors K'^-'^ = (Ak^-') , q^-'^) that also span the entire D- 
dimensional vector space. One can then find the dual 
basis A'^-') — (a(^',b(^^), expanded into 2-dimensional 
tiling vectors a^-'^ and {D — 2)-dimensional extensions 
b^-'^ using a D-dimensional orthogonality condition 

A« •K(J') = 27r*y, (4) 

which implies a Z?-dimensional completeness relation 

D 

^A(f")i^(^")=27r<5^., (5) 

thus satisfying the Ho condition Q . Note that the choice 
of the D—2 vectors {qj^^ , . . . , qjf^ ) is not unique, but once 
these are chosen the tiling vectors a^-') and their {D — 2)- 
dimensional extensions b'^-'^ are uniquely determined by 
the orthogonality condition Q). 

It can be shown |^ Eq. 5.5] that if p{r) is a sum 
of delta functions, at the positions determined by the 
dual grid method, then its Fourier transform p(k) is 
non-vanishing at most on the lattice of integral linear 
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combinations of the Ak^-'). If the Ak*^-'^ are hnearly 
independent over the integers this immediately gives 
us Bragg peaks at the positions required for multiple 
phase-matching. If the mismatch vectors are integrally- 
dependent some of the Bragg peaks may become extin- 
guished. To avoid this situation one may simply choose 
an independent subset of mismatch vectors to generate 
the structure, whose Fourier transform will then contain 
all the required Bragg peaks. Alternatively, one may re- 
fine the procedure for selecting the tiling vectors. We 
shall consider the latter option in greater detail else- 
where. The Fourier coefficient at k is given by 

p(k) = iVe"if / dse^'i-, (6) 



w 



where the sum — required only when the mismatch vec- 
tors are linearly-dependent over the integers — is over all 
q = X] "^i such that ^mj Ak^-'^ — k, the shift vec- 
tor f = X] /jb*--'^ , W is the so-called {D — 2)-dimensional 
window, given by the set of all points Ajb'^-'-' with 
< Aj < 1, s is a (-D — 2)-dimensional integration vec- 
tor, and V is the volume of the primitive cell of the D- 
dimensional real-space lattice generated by the vectors 

Methods for modulating the nonlinear coefhcient have 
been developed for many x^^^ materials, andjn particular 
for ferroelectrics and semiconductors . Although 
typically the nonlinear coefhcient is modulated only in 
one direction, the modulation methods are based on pla- 
nar technologies, hence 2-dimensional modulation is pos- 
sible 0- In ferroelectrics such as LiNbOs and KTiOP04, 
the modulation is based on reversing the electrical do- 
mains in the z direction of the crystal, and thus provides 

(2) . 

either a positive or negative value of the X33 coefficient. 
In the following calculation, we have assumed a mini- 
mum domain size of 2 /im, which can be achieved using 
the electric field poling technique Thus the actual 

nonlinear crystal g{v) is extended over a finite area 
and consists of positive domains of a given shape S, typ- 
ically small circles 8] or polygons such as hexagons [3 
or squares JOj. positioned at the vertices of the tiling in 
a negative background. Using the convolution theorem, 
its actual Fourier transform .g(k) is then 



5(k) = Ax (^p(k) ® J e'^^ 'd^rj J e^^ 'd'r, (7) 

where Ax is the absolute difference between the positive 

(2) 

and negative values used for X33 , p(k) is the sum of delta 
functions whose amplitudes are given in ©, and is 
the convolution operator. For example, if S* is a circle of 
radius R and A a rectangle of sides x Ly, the Fourier 
transform is 

Ji{kR) 



5(k) = gik,,ky) = 2SAAx 



p(k) 



smc 



kR 

kxLx 



where k = |k|, and Ji is the first Bessel function. There 
is, of course, much room for fine tuning of the design 
process — especially in the choice of the domain shape S — 
but we have clearly demonstrated that the required mis- 
match wave vectors indeed appear in the Fourier trans- 
form of the nonlinear photonic quasicrystal. Note in par- 
ticular for the circle, that care must be taken in the choice 
of the radius R such that the first zero of the Bessel func- 
tion, occuring at kR ~ 3.8317, does not extinguish the 
required Bragg peaks. 

As an example for our design procedure consider the 
color fan, depicted in Fig.^ whose input is a fundamental 
beam at frequency u> in the i, direction and whose out- 
put consists of the second, third, and fourth harmonics 
in different directions. Note that the pair of interacting 
beams, in all three frequency conversion processes, u and 
2ci;, are travelling in the same (x) direction to allow for 
efficient interaction throughout the full length of the pho- 
tonic crystal. Using the experimental values for the dis- 
persion in lithium niobate at 150°C, given in Ref. Eil and 
taking an input beam of wavelength A = 1550nm, the 
required mismatch vectors and the corresponding tiling 
vectors, satisfying the Ho condition Q, are given in Ta- 
ble m In order to keep all three mismatch vectors to 
be of the same magnitude we use for the first process 
{lj -\- Lo ^ 2lij) a third order quasi-phase-matching condi- 
tion, as described for example in Ref . 0. A section of the 
resulting photonic quasicrystal, produced with these pa- 
rameters using the dual grid method, is shown in Fig. |21 
Its Fourier transform is shown in Fig. O clearly show- 
ing Bragg peaks exactly at the required positions of the 
mismatch vectors. 
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FIG. 2: Structure of the photonic quasicrystal, produced with 
the parameters of Table U Scale is in ^.m. 



The method described in this Letter enables to achieve 
efficient frequency mixing of any arbitrary set of in- 
put waves, in order to generate a chosen set of output 
waves with specified wavelengths and propagation direc- 
tions. We have shown the implementation of the dual- 
grid method for a two-dimensional modulation of the 
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Mismatch vector 


Magnitude in (iim) ^ 


Direction in radians 


ff(Ak(^)) 


Ak'i^ = Akii^a 


1.0183 





0.0404 


Ak(2) = Ak22^4 


2.8179 


0.0456 


0.0431 


Ak(=" = Aki2^3 


0.9530 


-0.0966 


0.1671 


Tiling vector 


Magnitude in /im 


Direction in radians 






10.330 


-1.4801 






18.234 


1.5001 






42.297 


-1.5289 





TABLE I: Calculated mismatch vectors for the color fan and the corresponding tiling vectors. The right-hand column lists the 
analytically-calculated Fourier coefficients for the three processes assuming an infinite structure and circular positive domains 
of radius 2.9212^m in a negative background. The corresponding values obtained in a 1-dimensional periodic crystal for first 
and third order processes are 2/-K ~ 0.6366 and 2/37r ~ 0.2122, respectively (3, p. 2635]. 




FIG. 3: Numerical Fourier transform of the crystal in Fig. |21 
showing Bragg peaks at the positions of the required mis- 
match vectors, as indicated by arrows. Scale is in {fim)~^ . 

nonlinear coefScient. The method is, of course, appH- 
cable in any number of dimensions. The implementa- 
tion in one dimension is straightforward and requires that 
the wave-vectors of the interacting waves, and thus the 
mismatch wave-vectors, are all collinear. For example, 
the one-dimensional Fibonacci structure ^| can easily 
be generated by choosing two collinear mismatch wave- 
vectors having a ratio of r = (1 -f V5)/2 between them. 

An interesting device that is enabled by the proposed 
method is a nonlinear mixer that can generate the same 
harmonic at several different directions. This mixer can 
be highly useful for cavity-enhanced nonlinear processes, 
such as resonant second harmonic generation and opti- 
cal parametric oscillation, with the unique feature that 
all propagation directions in the cavity contribute to the 
nonlinear process. In addition to the possibility to gener- 
ate multiple harmonics at multiple directions, as shown 
here, some interesting applications of this method can 
benefit from the possibility of generating large nonlin- 
ear phase shifts by cascaded nonlinear processes. These 
applications include generation of multicolor optical soli- 
tons, as well as all-optical deflection and switching of 
light |10|. 



This research is supported by the Israel Science Foun- 
dation through Grant No. 278/00, and by the Israeli Min- 
istry of Science. 



* Electronic address: lronlif@tau.ac.ill 
[1] J. A. Armstrong, N. Bloembergen, J. Ducuing, and P. S. 

Pershan, Phys. Rev. 127, 1918 (1962). 
[2] J. A. Giordmaine, Phys. Rev. Lett. 8, 19 (1962). 
[3] P. D. Maker, R. W. Terhune, M. Nisenoff, and C. M. 

Savage, Phys. Rev. Lett. 8, 21 (1962). 
[4] N. Bloembergen and A. J. Sievers, Appl. Phys. Lett. 17, 

483 (1970). 

[5] M. M. Fejer, G. A. Magel, D. H. Jundt, and R. L. Byer, 
IEEE J. Quantum Electron. 28, 2631 (1992). 

[6] V. Berger, Phys. Rev. Lett. 81, 4136 (1998). 

[7] N. G. R. Broderick, G. W. Ross, H. L. Offerhaus, D. J. 
Richardson, and D. C. Hanna, Phys. Rev. Lett. 84, 4345 

(2000) . 

[8] S. Saltiel and Y. S. Kivshar, Opt. Lett. 25, 1204 (2000). 
[9] X.-H. Wang and B.-Y. Gu, Euro. Phys. J. B 24, 323 

(2001) . 

[10] S. M. Russell, P. E. Powers, M. J. Missey, and K. L. 

Schepler, IEEE J. Quantum Elect. 37, 877 (2001). 
[11] S.-N. Zhu, Y.-Y. Zhu, and N.-B. Ming, Science 278, 843 

(1997). 

[12] S.-N. Zhu, Y.-Y. Zhu, Y.-Q. Qin, H.-F. Wang, C.-Z. Ge, 
and N.-B. Ming, Phys. Rev. Lett. 78, 2752 (1997). 

[13] K. Fradkin-Kashi and A. Arie, IEEE J. Quantum Elec- 
tron. 35, 1649 (1999). 

[14] K. Fradkin-Kashi, A. Arie, P. Urenski, and G. Rosenman, 
Phys. Rev. Lett. 88, 023903 (2002). 

[15] Y. S. Kivshar, A. A. Sukhorukov, and S. M. Saltiel, Phys. 
Rev. E 60, R5056 (1999). 

[16] S. M. Sahiel, A. A. Sukhorukov, and Y. S. Kivshar, 
Progress in Optics 44 (2004), (nlin.PS/0311013). 

[17] M. Senechal, Quasicrystals and geometry (Cambridge 
University Press, Cambridge, 1995). 

[18] R. Lifshitz, Found. Phys. 33, 1703 (2003). 

[19] N. G. de Bruijn, Proc. K. Ned. Akad. Wet. Ser. A 84, 39 
(1981). 

[20] F. Gahler and J. Rhyner, J. Phys. A 19, 267 (1986). 
[21] D. A. Rabson, T.-L. Ho, and N. D. Mermin, Acta Crys- 
tallog. A 44, 678 (1988). 



5 



[22] D. A. Rabson, T.-L. Ho, and N. D. Mermin, Acta Crys- 

tallog. A 45, 538 (1989). 
[23] T.-L. Ho, Phys. Rev. Lett. 56, 468 (1986). 
[24] M. Yamada, N. Nada, M. Saitoh, and K. Watanabe, 

Appl. Phys. Lett. 62, 435 (1993). 



[25] L. Eycrs, P. ,L Tourrcau, T. ,J. Pinguet, C. B. Ebcrt. 
J. S. Harris, M. M. Fejcr, L. Becouarn, B. Gerard, and 
E. Lallicr, Appl. Phys. Lett. 79, 904 (2001). 

[26] D. H. Jundt, Optics Lett. 22, 1553 (1997). 



